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Abstract 

We present the calculation of the inclusive Xi^-distributions of charmed hadrons, 
produced in high-energy E~-beam. The calculation is based on the modified 
mechanism of charmed quarks fragmentation as well as on the mechanism 
of c-quark recombination with the valence quarks from initial hadrons. We 
predict the additional asymmetry in the production of charmed hadrons due 
to the different distributions of the valence s and d quarks in S~-beam. 
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Introduction 



In the hadronic production of particles with open charm an interesting sit- 
uation, connected with the interaction of charmed quarks with the hadronic 
remnant, takes place (see [|l], ||, || and references therein). 

Let us remind that such a problem does not arise in the case of e^e~ 
annihilation, where a heavy quark is hadronized due to its own radiation and 
all the description is reduced to the fragmentation function: 

= D{z, Q^) (1) 

here z is the part of the c-quark momentum, taken by charmed hadron and 

is the square of total energy in the e~^e~ annihilation. 

It is well known that a simple transition of factorization approach in the 
form of (|l|) to the charmed particle production in hadronic interaction leads to 
the valuable deviation from experimental data. Indeed, there is a substantial 
difference in the yield of different charmed hadrons in the fragmentation 
region — >■ 1, x = 2pH/y/s) of initial hadrons. 

This charmed particle yields asymmetry A (or "leading particle ef- 
fect") is defined as follows: 

^ ^(/ea(im(7) — ^{nonleading) 
^{leading) + ^{nonleading) 

Here symbol 'leading' or 'non-leading' refers to charmed hadron with or with- 
out light quarks identical to valence quarks from initial hadrons. For example, 
in 7r~(ud)-beam the mesons D^{cu) and D~{cd) are leading particles, while 
Z)°(cu) and D~^{cd) should be considered as non- leading ones. 

There are many theoretical articles devoted to description of a such phe- 
nomenon IQ. In our previous 0, ^ publications we have taken into account 
the interaction of c-quarks with quarks of the initial hadrons and obtained a 
good agreement with the data in the spectra description at large x and the 
x-asymmetry dependence in n'N and S~A^-interactions. 

The recent data on the charm production in the E~-beams [0, ^ are 
of special interest due to the fact that the beam hadron (here S) has s- 
quark and as consequence the distribution functions of the valence quarks 
in the S~-baryon differ from that of proton U^. It has to lead to the set of 
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observable effects and, in particular, to the another asymmetry in the yield 
of the charmed hadrons that in the case of proton-proton collisions. 

Preliminary data presented by the SELEX Collaboration on the charm 
yield in the S~p and 7r~p collisions show however some disagreement with 
other experiments and general theoretical assumptions, in particular, in 
the expected charmed particle and anti-particle yields in the central region. 

In this article we consider very briefly this problem from general the- 
oretical point of view (Section 1). In the Section 2 we give a short de- 
scription of the fragmentation mechanism. Our approach for calculation of 
x-distributions of charmed-particle production in recombination mechanism 
is given in the Section 3. We present the predictions for differential distribu- 
tions of charmed hadrons produced in E~A^-interactions in the Section 4. A 
short summary of the results is given in the Conclusion. 



1 Asymptotic behavior in the central region 

Before discussing in detail of the spectra behavior and the charge asymme- 
try in the hadron yield let us consider general behavior of charmed particle 
spectra. 

According to the generalized optical theorem ||^ an inclusive spectrum 
is connected with the discontinuity of the sixth particle amplitude by the 
following relation (see Fig. 1.): 

where s is the total energy squared. 

Let us introduce two invariants (see Fig.l) 

At the high energies the asymptotic behavior of the Ms^s amplitude as a 
function of these invariants si,2 is determined by leading Regge trajectories, 
namely, the Pomeron {V) and secondary trajectories (T^j), related to p, cj, 
/, A2 mesons. In so doing, the VV, TZV, and TZTZ contributions become 
important. 
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In the central region (x ~ 0), where the kinematic invariants si and S2 
are large, one has 

si ^ ^/sm±e~'^ , S2 ^ y/sm^e"^ , (4) 



where y is the rapidity and m± = y + p j . 

Double Regge representation is a good approximation for the Ms^s amphtude 
(see Fig.l) in this kinematic region. So, one has 

E^-^Y.ksTs? (5) 



where and are the intercepts of the leading Regge-trajectories |g, [11 
For Pomeron trajectory one has ap ~ 1, while the f,p,uj,A2 trajectories 
have the intercept aR ~ 1/2. 

The fij{m±) are unknown functions of m±, but it seems that these func- 
tions are universal and do not depend on the type of produced particles ||, |]. 
These functions can be determined by fitting to experimental data. All the 
dependence on quantum numbers of initial and final particles is determined 
by the coupling constants of the secondary Reggeons 71, included to the 
definition of the fij{m±) functions. 

After substituting Si^2 from (|) the expression takes the form: 

This equation provides a good description of the transition to the asymptotic 
regime in the "central" region for all variety of the yields of tt- and K- 
mesons 0. 

For the D- and S-meson cases the expression (^) simplifies due to the 
fact that the contribution of the trajectories, connected with a charmed c- 
quark in D-meson (J/ ip, Xc trajectories) or beauty 6-quark in i?-meson (T, Xb 
trajectories), is strongly suppressed by the Zweig rule for the coupling with 
initial mesons or nucleons. Therefore, in the sum of Eq. (j^) there are no 
contributions with two secondary Reggeons (7?.7^-contributions) in the upper 
and bottom "shoulder" of the diagram in Fig. 1. So, one gets more simple 
expression: 

e'^^ ^ fvvirn^) - j^{fnvim^)e^/' + hnMe"^"}. (7) 
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Note, that the contribution of the first VV-ieim. in the equation (|^ is 
the same for the particle and anti-particle and is proportional to the total 
cross-section of the interaction of initial hadrons. Thus, the normalized cross- 
section 

1 ^d?(y{hh ^ DX) 

at high energy limit has one universal limiting value independently on the 
type of colliding hadrons. 

For the case of the Dg-meson production in the S or i^-beams the leading 
intercept is related to 0-meson, aj = a,/, ~ and there is no contribution 
from the third PT^.-term in Eq. (|^. The 0-trajectory can be connected with 
S-particle only, because there are no valence strange quarks in the proton. 
As a result, one has only two terms: 

— ^^^^ ^ fM^±) - ^kv{m^)e' (8) 

The contribution of secondary trajectories determines the difference be- 
tween the yields of the particle and anti-particle in the central region. For 
instance, for pp and pp-collisions one has 



a y b y 



where the a and b coefficients depend on both combinations of secondary 
Reggeon couplings and particle type of D- or 5-meson. 

In the most experiments on fixed target the point y = (x = 0), as a rule, 
is out of reach due to the experiment conditions, and the measured spectrum 
starts from x ~ 0.1. With a simple substitution y ~ ln(a:^) for D-mesons 
the equation takes the form 

E—r. frp{m±) - {fnvijn^)^/ml^x + /p7e(mx)-^^^} . (9) 

a tot a V y sx 

For Dg-meson production in the S or beams one has: 

—E ^ o-ps) ^ f^^(rn±) - fTivim±)m±x (10) 
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The asymmetry in the particle yield, which depends on the quantum num- 
bers of initial hadrons and observed charmed particles, is determined by the 
second term in these expressions, which is different for the D and Dg cases. 
One can see that the transition to the asymptotics at fixed x is achieved 
much faster than in the regime of the fixed rapidity. It is seen also that the 
behavior in the vicinity of a; = is determined by the value of intercept of 
secondary trajectory connected with valence quark, common for the beam 
hadron and observed particle. 

It is quite evident that the applicability region of the above expressions 
is severely restricted by the necessity to fulfill the condition of large values of 
the Si and S2 invariants in the case of the Regge approximation. However, the 
general conclusions on the character of an asymptotic behavior of charmed 
particle spectra, which can be obtained from them, on the one hand, are 
quite definite, and, on the other hand, agree well with the parton model 
predictions. 

We wish to stress once more, that from general theoretical consideration 
one should expect equal yields of the charmed particle and anti-particle in 
the central region [x ~ 0) when s —>■ 00: 

E^^U^o-E^^l^^o (11) 
As we will show such behavior is agree also with parton model predictions. 



2 Fragmentation mechanism 

In our previous publications ^ we developed the phenomenological model 
model, where the hadronization of charmed c-quark is described by the sum 
of two mechanisms, namely 

dx dx dx ' 

where the first term corresponds to the charmed quark fragmentation, while 
the second term takes into account the charmed c-quark interaction with 
valence quarks from initial hadrons (recombination). 
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In the fragmentation mechanism the inclusive cross section for the charmed 
hadrons (D-mesons) production has the form as follows: 

„ d^a^ f ^ d^a(hih2 cX) ^, , ^ s , . 

""1^ = J 4v — ' - (1^' 

where z = \ph\/\Pc\ is the fraction of the c-quark momentum carried away 
by the charmed hadron H. 

The parameterization of the fragmentation function D{z) (for example, 



in the form from [|T2[ or [0) can be found by fitting to the data from reac- 
tion e~^e~ D{cq) X. However, in hadronic collisions the situation is more 
complicated. Indeed, the use of the fragmentation function is justified at 
asymptotically large values of the invariant mass of cc pair or high pj. How- 
ever, in hadronic production of the charmed particles the main contribution 
into inclusive charm production cross section results from c-quarks with low 
values of the invariant mass of cc pair [M^c > 2mc) |10[ . These quarks domi- 
nate in the region small x. At the same time there is large amount of partons 
from initial hadrons in the same (central) region of x. Therefore, the c-quark 
in combination with one of such a parton can easily to produce a charmed 
hadron. Such a process occurs practically without any loss of c-quark mo- 
mentum (i.e. ph ~ Pc)- Therefore, in the small x region one should expect 
the coincidence of the spectra of charmed hadrons and c-quarks. Whereas at 
high X region one may use the conventional fragmentation mechanism. 

Following these arguments we have proposed the modified form of the 
fragmentation function 

D''^{z,M,,)r^z-"^^^-^\l-z), (14) 

with two additional conditions on a(Mcc): 

a{Mcc) -oo D{z) (5(1 - z) at M^c 2mc, 



a{Mcc) -2.2 L)(z) ^ ^-"-(1 - z) at M^c ^ Mq. ^^^"^ 

The explicit form of a{Mcc) equal 

1 - 3/i(M,5) 



1 - n{Mcc) 



where 



In go 
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Note, that the usage of the fragmentation function assumes the absence 
of the interaction of the produced heavy c-quark with the remnants of the 
initial hadrons. Therefore, it should be no difference between the spectra 
of charmed and anti-charmed hadrons. Moreover, any modification of the 
fragmentation mechanism can not reproduce the production asymmetry (the 
leading particle effect). 



3 Recombination mechanism 

The fragmentation mechanism can be apply for the production of the cc pair 
in the color-singlet state or for high production of the open charm. On 
the other hand, for the case of the hadronic production of color cc pair with 
small pt one should takes into account the possibility of charmed c and c 
quarks interaction with the initial hadron remnants. Therefore, due to the 
different valence quarks in the initial hadrons one may expect the different 
inclusive spectra of the final charmed hadrons. 

In the parton model framework, a heavy c-quark should interact with a 
high probability with its nearest neighbor in the rapidity space able to form a 
color-singlet state with it. In some cases, the heavy anti-quark (quark) may 
find itself close (in rapidity space) to a valence light quark (diquark) from 
the initial hadron. This would result in the formation of a fast heavy meson 
(baryon) in the fragmentation region of the initial hadron. 

We use the model |^ to describe the production asymmetry for charmed 
hadrons. In this model the interaction of the charmed quarks with valence 
quarks from the initial hadrons describes with the help of the recombination 
function p, ^. The recombination of the valence qv and c quarks into D- 
meson is described by the function of -Rm( 



where C,q = Xq/x and = z/x, while Xq, z, and x are the fractions of the 
initial-hadron cm. momentum that are carried away by the valence g-quark, 
charmed c-quark, and the D(cg)-meson, respectively. The corresponding 
recombination of three quarks into baryon can be described by means of the 
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similar recombination function: 



Rb{xi,X2,z;x) 



r(3 — ai — a2 — Oc) 



(17) 



r(i-ai)r(i -a2)r(l-a,) 



These functions take into account the momentum conservation and the prox- 
imity of partons in the rapidity space. Actually, the recombination function is 
the modulus squared of the heavy meson (baryon) wave function in momen- 
tum space, being considered in the infinite momentum frame in the valence 
quark approximation. 

With the help of the function R{xv, z; x) describing the recombination of 
the quarks qv and c into a meson, we represent the corresponding contribu- 
tion to the inclusive spectrum of D meson as follows 0: 



where x* = lEj^ and x = 2pi/y/s (here, E and pi are the energy and 
longitudinal momentum of the D-meson in the c.m.s. of the initial hadrons); 
xy and z are the momentum fractions carried away by the valence quark and 
heavy anti-quark, respectively, and xyz* ^f^"^^ is the double-differential cross 
section for the simultaneous production of the quarks qv and c in a hadronic 
collision. The equation describing the production of charmed baryon has an 
analogous form. 

The parameter Rq is the constant term of the model, that determines the 
relative contribution of recombination. We fit the data on charmed hadrons 
production in n~N collisions and found that Rq ~ 0.8 0. 

Note, that using of the recombination with the valence quarks provides a 
well description of the leading particle effect. At the same time its contribu- 
tion in the total inclusive cross section production of the charmed particles is 
sufficiently small (~ 10%). This mechanism dominates in the high x region. 

Note, that this model provides more or less successful description of the 
charmed D-meson production in vr" A^-interactions (see Fig. 2, 3 and p for 
details). 




(18) 
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4 Charm production in S— beam 



In this section we consider the charmed hadron production in high-energy 
beam of S~-hyperon. First of all we expect the various behavior of the 
distributions of the valence d- and s-quarks. Indeed, as a first approximation, 
the distribution of valence quark in the baryon -B(gig2'?3) can be presented 
as follows 

V^^{x) oc - (19) 

where is the intercept of the leading Regge-trajectory for g^-quark, while 
7b ~ 4. Note, that due to violation flavor S'f/(A^)-symmetry, we have differ- 
ent intercepts for d{u)- and s-quarks [|ll|, [T^: 

au = ad = ^, as ^0, a^^ -2.2. (20) 

As a result, the x-dependence of the valence d and s-quark in the H'^sdd)- 
hyperon has the form as follows 0: 

vf~^(l-x)3•^ v^'^ii-xf (21) 

V X 

It is seen from the ( |2ll) that the valence s-quark in the S~-hyperon has 
harder x-distribution than that for d-quark. Note, that gluon distribution 
in S-hyperon has the same form as in usual nucleon. As a result, we expect 
more harder x-dependence of s-quark distribution in S-hyperon then u-quark 
distribution in nucleon. Thus, one may expect slightly harder spectrum of 
D^-mesons, produced in S-beam, then spectrum of D°-mesons, produced in 
nucleon-heam. On the other hand, d-quark in E-hyperon should be slightly 
softer then that in nucleon. As a result we should observe the different 
Xi?- dependence of spectra of charmed hadrons with d- or s-quarks, namely, 
D^icd) and D^{cs), S°(c(is) and S^(c(i(i), etc. 

We use LO formulas for calculations of the cross sections for quark- 
antiquark and gluon-gluon annihilation into cc-pair. We set = 1.25 GeV 
and the value of strong couphng constant equals 0.3. Then for the cross 
section value of S~ p interaction at Plab = 600 GeV one has: 

(t(S-j9 ^ ccX) lib (22) 

In our calculations we do not pretend to reproduce the absolute value of 
this cross section (see [|1^], for detail consideration of this problem). We 
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concentrate on the description of Xi^-distribution of charmed mesons and 
baryons. 

Recently, we have done the calculations of asymmetry of x-spectra for 
Dg and Ac hadrons produced in S~ beam with the energy of 340 GeV. 
The WA89 Collaboration have compared these predictions with their exper- 
imental data 0. This comparison is presented in Fig. 4. As is seen from the 
picture the asymmetry in D-meson production in S-beam is different from 
that in 7r~-beam and there is more pronounced asymmetry for D^-meson 
production. 

Below we present the predictions for charmed hadrons production in E~- 
beam at 600 GeV energy. The corresponding distributions (integrated over 
Pt) are presented in Fig.5. We may see from these figures, that the considered 
charmed quark interaction in the final state (recombination) leads, indeed, to 
noticeable differences in XiT'-spectra. These differences can be explicitly seen 
in Fig. 6, where we present the corresponding asymmetry A (see (0) for defi- 
nition). The most non-trivial prediction of the proposed model is presented 
in two lower plots in Fig.6, where we present the ratio of the inclusive spectra 
of Dl~{cs) and D*~{cd) mesons as well as and S° baryons. Indeed, due 
to the difference of the valence d and s quarks in the S~-beam (see (pT])) we 
expect the additional asymmetry in the leading charmed particles produc- 
tion. 



5 CONCLUSION 

In the present article we wish to stress once more, that the source of the 
observed asymmetry in charmed hadron production is the interaction of pro- 
duced charmed quarks with valence quarks from initial hadrons. Note, the 
model under consideration provides also the additional method to measure 
the valence quark distribution functions of S-baryons. 
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P.S. When this article was finished we received the article with the 
recent results of E791 Collaboration. They present the inclusive x and pt 
distributions of charmed hadrons, produced in 7r~-beam. In particular, they 
observed a noticeable asymmetry in and A~ baryons production in the 
forward region (about 13%). Note, that our model can not explain this result. 
Due to the equal amounts of valence u and d quarks in the 7r~-beam we 
expect an equal yields of the A+(cMd) and A.~{c\id) baryons in the 7r~-beam. 
Moreover, the most theoretical models predicts also zero asymmetry for 
these particle production. 
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Figure 1: Six particle Ma^s amplitude. 
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Figure 2: Differential distributions ^ for the energy of E'tt = 
250 GeV The dotted (dashed) histogram corresponds to the re- 
combination (fragmentation) contribution. The sohd histogram rep- 
resents their sum. The cross sections are presented in fih (see Q for 
details). 
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Figure 3: The description of the asymmetry A(x) in tt p collisions 
(see |6| for details) j^g 




Figure 4: The asymmetry in the charmed hadrons produced in S p 
interactions at Plab = 340 GeV g. 
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Z p PL.B = 600GeV 
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Figure 5: ^-distributions of cliarmcd mesons {D*^ and and 
baryons produced in S~ p interactions at Plab = 600 GeV. The solid 
(dashed) curves correspond to leading (non-leading) charmed hadron 
production. 
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r p PjAB = 600 GeV 




Figure 6: Production asymmetry A{x) from @ for D*"^ and D*^ 
mesons (two upper histograms), produced in T,~ p interactions at 
Plab = 600 GeV. Two lower plots present the ratio of the spectra of 
charmed hadrons with and without strange quarks (i.e. the D*^ /D*^ 
and S[?(cds)/S[!(cd(i) ratios). 
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